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Abstract 

Non-Abelian geometric phases can be generated and detected in certain 
superconducting nano circuits. Here we consider an example where the 
holonomies are related to the adiabatic charge dynamics of the Josephson 
network. We demonstrate that such a device can be applied both for adia- 
batic charge pumping and as an implementation of a quantum computer. 



* permanent address 



1 



If a quantum system is prepared in a superposition of two states, a physical observable 
associated with this system can exhibit oscillatory behavior depending on the relative phase 
of the two states. Interference can be induced during the dynamical evolution of the system, 
in this case we refer to the accumulated phase as the dynamical phase. Interference can also 
be of geometrical nature if the parameters (coupling constants, external fields,. . .) of the 
Hamiltonian are varied cyclically After Berry's original work, considerable attention has 
been devoted to the interpretation, generalization, and detection of geometric phases ||. 
An important generalization is when the adiabatic cyclic evolution involves a degenerate 
eigenspace of the Hamiltonian. In this case it has been shown by Wilczek and Zee Q 
that the evolution over a closed path does not result in a phase change but it leads to a 
superposition of the degenerate eigenstates and the geometric phase acquires a non-Abelian 
structure. Originally investigated in Nuclear Quadrupole Resonance ||, more recently it 
was shown that non-Abelian phases occur in the manipulation of trapped ions |5||J. 

Apart from its fundamental importance, geometric interference has interesting applica- 
tions in the field of quantum information processing Implementations of quantum 
computers so far include optical systems and liquid-state NMR |Sj as well as solid-state 
devices based on superconductors [[HJ and on semiconductors |TI|] . Recently it has been 



shown that quantum computation can also be implemented by geometric means (geometric 



quantum computation) using Abelian [12fl as well as non-Abelian [T3||| phases. 

Non-Abelian phases can also appear in the quantum dynamics of superconducting 
nanocircuits |14| |, this is what we will show in this work. There are various interesting 



aspects associated with this analysis. In addition to their possible detection, which is in- 
triguing by itself, the existence of non-Abelian phases in superconducting nanocircuits leads 
to a new scheme for adiabatic charge pumping and allows to implement solid state holo- 
nomic quantum computation. Some parts of our proposal are on purpose speculative. The 
adiabatic manipulation of degenerate subspaces and the degeneracy condition itself is non- 
trivial to achieve for an artificially fabricated device. We believe, however, due to the rapid 



development in the control of artificial two- level systems in solid-state devices JL5| , the real 



ization of geometric interference in mesoscopic systems has become plausible. Furthermore 
the applications of non-Abelian phases in pumping and computation are interesting new 
directions to pursue. 

In our discussion of non-Abelian phases in Josephson junction circuits we follow the spirit 
of the schemes described in Refs. The starting point is the network shown in Fig.|TJa) . It 
consists of three superconducting islands labeled by j = L,M,R (Left, Middle, Right) each 
of which is connected to a fourth (Upper) island labeled with U. Gate voltages are applied to 
the three bottom islands via gate capacitances. The device operates in the charging regime, 
that is the Josephson energies Jj (j = L, M, R) of the junctions are much smaller than the 
charging energy E c of the setup. Each coupling is designed as a Josephson interferometer (a 
loop interrupted by two junctions and pierced by a magnetic field) as shown in Fig.|l]a. Thus 
the effective Josephson energies Jj can be tuned by changing the flux in the corresponding 
loop. Electrostatic energies can be varied by changing the gate voltages V g . 

Let us first analyze the electrostatic problem (i.e. Jj = 0). For the sake of sim- 
plicity we assume that all capacitances are equal to C and we consider identical gate 
charges q g = C g V g /{2e) for the three bottom islands. The charge states are indicated as 
| njj, ni, n M-, tir) where n labels the number of Cooper pairs in the corresponding island | |16| . 
For gate charges q g ~ 1/2 and 1 < 2qu + Sq g < 2 (where qv = CjjVu/{2e), see Fig.|l|a), only 
four charge states are important as long as T Ec = (2e) 2 /(4C). Three of these charge 
states, | 0,1, 0,0), | 0,0, 1,0), | 0,0,0,1) are degenerate. Their charge configuration corre- 
sponds to one excess Cooper pair in one of the islands j = L, C, R, and none in the island 
U. The fourth state | 1, 0, 0, 0) has one excess pair on the island U and none on the other 
islands. All other charge states are much higher in energy. 

The Josephson couplings Jj allow for tunneling between the upper island and each of 
the bottom islands. The quantum-mechanical Hamiltonian of this simple four-state system 
reads (in complete analogy with Refs. 0f|) 
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H = 5E C \1, 0, 0, 0) (1, 0, 0, 0| + i[J L ($ L )|l, 0, 0, 0) (0, 1, 0, 0| + 

+J M ($ M )|1,0,0,0)(0,0,1,0| + J fi ($«)|l,0,0,0)(0,0,0,l| + H.c] (1) 
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where 5Ec = —Ec 
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is the energy difference between the three degenerate 
states and the fourth one, and $j = <&j/$o are t ne external magnetic fluxes in units of the 



flux quantum $o = hc/2e [ 17fl . In all the manipulations described below the gate voltages 
will be kept fixed. The three fluxes {5>l, $m> ^/j} are the parameters which will be varied 
cyclically. In general all the SQUID loops could be asymmetric, although it is not necessary 
for the purpose of our discussion. 

The Hamiltonian defined in Eq.(|l]) can easily be diagonalized. The lowest and highest 
eigenstate are non-degenerate. The peculiar feature, exploited in is that the other 
two states (with zero energy) are degenerate for arbitrary values of the couplings Jj. The 
subspace is spanned by the eigenstates (not normalized) 

|£>i> = -Jm|0, 1,0,0) +J L \0, 0,1,0) , 

\d 2 ) = -j r (r L \o, 1,0,0) + ^10,0, i,o)) + (|j L | 2 + 1 j Af | 2 )|o, 0,0,1). (2) 

By manipulation of the external magnetic fluxes it is possible to generate non-Abelian phases. 
We will show that by means of such phases adiabatic charge pumping and holonomic quan- 
tum computation can be realized with superconducting nanocircuits. 

Charge pumping - For this purpose it is sufficient to have only symmetric SQUID loops. 
In contrast to the well-known turnstiles for single electrons or Cooper pairs pT8|JT9|j20[1 (in 
which the gate potentials are modulated periodically), here charge is transported through the 
chain (from the L-island to the R-island) by means of modulating the Josephson couplings 
while keeping the gate voltages unchanged. The pumping cycle goes as follows. The system 
is initially prepared in the |0, 1,0,0) state (i.e. the state \Di) with Jl,Jm = ) where the 
Cooper pair is in the left island. This can be achieved by turning off all Josephson couplings 
and coupling the L-island to a lead which provides the extra Cooper pair. Once the charge 
is on the island, one should change adiabatically the magnetic fluxes along a closed loop 7. 



At the end of the loop the initially prepared stated \Dx) will be mapped into the following 
rotated state: \Di) — > U 7 \Di), where the unitary matrix U 1 may be expressed as ||: 

d_ 

'~<j=L,M,R d ®3 

Path ordering P is required as, in general, the matrices Aj do not commute along the path. 
If the path 7 is chosen in the ($£, $r)- plane (at fixed 5>Af = 0), as shown in Fig.|l]b it can 
be shown that, after one adiabatic cycle, the final state of the system is |0,0,0, 1), i.e. one 



U, = Pex V <f £ Ajd$j ; {Aj) a>p = (D a \~\D p ), «, /? = 1,2. (3) 



Cooper pair has been transported through the chain of the three islands pi 

The mechanism described here relies entirely on the geometric phase accumulated during 
the cycle and can be generalized to describe pumping of a single Cooper pair through iV 
superconducting islands. The connection between pumping and geometric phases has been 
discussed by Pekola et al. ||20|| . The crucial difference is that here only the Josephson 
couplings have to be varied. During the cycle, exactly one Cooper pair is transported, in 
this sense there are no errors due to the spread of the wave function discussed in |20] . 



There are drawbacks though, mostly related to the fact that the degenerate states are not 
the ground state and relaxation processes may become important. 

Quantum computation, one qubit - The pumping process illustrated so far is nothing 
but one of the key elements to construct a quantum computing scheme using non-Abelian 
phases. Proceeding along the lines of Ref. ||, we point out the necessary ingredients and 
the differences which arise in the case of the Josephson junction setup. The nanocircuit 
presented in Fig.[TJa constitutes the qubit. The logical states to encode information in this 
implementation are 

10)^=10,1,0,0) 
\l) t = 10,0,0,1) . 

The other two charge states (|1, 0,0,0) and |0, 0,1,0)) serve as auxiliary states. To show 
that the implementation is possible, it is sufficient to provide explicit representations for 
the gates U\ = e 121 ' 1 '"' 1 ' and U2 = e iS2<7y , describing rotations of the qubit state about the 



z axis and the y axis, respectively. In this case only one asymmetric SQUID (as shown in 
Fig.[I]) is required to implement the one-qubit operations. 

The gate U\ is a phase shift for the state \l)i while the state |0)^ remains decoupled, 
i.e., Jl = during the operation. In the initial state we have Jr = 0, so the eigen- 
states { l-Di) , \D 2 )} correspond to the logical states {|0)^, ll)^}. The control parameters 
evolve adiabatically along the closed loop C\ in the ($m, $/j)-plane from = 1/2 
to <&r = 1/2 (see Fig.|2|a). By using the formula for holonomies Eq.(^|) one can show that 
this cyclic evolution produces the gate U\ with the phase Si : 



r sin (27t$r) 
E x = (7i & d$ M d®R- = - 



■M,\ 2 



(4) 



(\M$r)\ 2 + \Jm(® 

where «S(Ci) denotes the surface enclosed by the loop G\ in M.. and o\ = 4tt 2 Jr(0) 2 (| Jmi\ 2 — 

\Jmt\ 2 )- 

Similarly we can consider a closed loop C 2 (see Fig.^o) in the (<&£, $#)-plane at fixed $m = 0, 
and let the control parameters $l and undergo a cyclic adiabatic evolution with starting 
and ending point 5>l = 5>_r = 1/2. This operation yields the gate U 2 with phase E 2 



r _ sin ( 7r$i ) sin (hQr ) 
E 2 = a 2 f d$ R d$ L * '- * ^ 



—U2 ( 5 ) 

,5(C2) (Jm(0) 2 + \Jr($r)\ 2 + \J l ($l)\ 2 ) 

where S(C 2 ) denotes the surface enclosed by the loop C 2 in A4, and a 2 = 4ii 2 \ Jr(0)\ 2 (\ Jmi\ 2 + 
\Jmt\ 2 ) where we have assumed Jl(0) = Jr(0). Obviously the pumping cycle discussed above 
is a special case of the gate U 2 with E 2 = 7r/2. 

Quantum computation, two qubits - It turns out that it is possible to implement a con- 
ditional phase shift U3 = e jE 3|n>«<n| coupling two qubits via Josephson junctions. These 
junctions should be realized as symmetric SQUID loops such that the coupling can be 
switched off. The capacitive coupling due to these SQUID loops can be neglected if the 
capacitances of the junctions are sufficiently small [f22]| . 

By setting 5Ec = (this was not necessary in the one-qubit case) and by coupling the 
qubits as shown in Fig.^ja) we obtain the Hamiltonian 
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rr 1 

-n2qubit — g 



J$\UU)(UM\ + J x |ll)/(^| + H.c. (6) 



where we introduced the notation \U) = |1, 0,0,0) and \M) = |0, 0,1,0) for the aux- 
iliary states which are coupled by the interaction between the qubits. The matrix 
element J x = Jx($ur,$ru) is given by J x = -{1/2)J UR ($ UR ) J ru {$ur) V where 
/i = 1/ 5Eq~ + 1/ $Ec . Here 5E^~ and SE^ + denote the charging energy difference be- 
tween the initial and the intermediate state (see below). The coupling is of second order in 
the Josephson energies since the inter-qubit coupling junctions change the total number of 
pairs on each one-bit setup. Thus the coupling occurs via intermediate charge states which lie 
outside the Hilbert space of the two-qubit system. These are states, e.g. |0, 0, 0, 0)®|1, 0, 0, 1), 
without excess Cooper pair on the first qubit and two excess pairs on the second qubit. We 
have abbreviated the charging energy difference between the corresponding state and the 
initial qubit state by SEq + , and we have denoted the external magnetic fluxes in the cou- 
pling SQUID loops by <&ur and $ru- 



While Jxi^uRi &ru) is the only off-diagonal coupling of second order, there are also second- 
order corrections of the diagonal elements, i.e. of the energies of the two-qubit states. These 
corrections would lift the degeneracy and thus would hamper the geometric operation which 
is based on the degeneracy of all states. It is therefore crucial that it is possible to compen- 
sate these corrections and to guarantee the degeneracy. It is easy to see that by adjusting the 
gate voltages the energy shifts can be canceled. Note that during the geometric operation 
the values of the Josephson couplings are changing and therefore also the energy shifts are 
not constant. Consequently their compensation by means of the gate voltages has to follow 
the evolution of the parameters. 

Let us now show explicitly how the gate U3 can be achieved. To this aim, we consider a 
closed loop C 3 in the (<3>y , <3>[/#)-plane at fixed <&ru = 0. (See Fig.|^b). If the control pa- 
rameters §ur and & M undergo a cyclic adiabatic evolution with starting and ending point 
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®ur = 1/2, $m = 0, ^ e geometric phase obtained with this loop is 

r -12) - sin ( 27T ®ur) 

£3 = a 3 f d&$d$ UR * 

Js{Cs) (l4M?)l 2 + I* 2 Jru(0) 2 \Jur(*ur) 

with a 3 = Att 2 J ru (0)\\J 2 mii - \J Mr \ 2 ) and J UR (0) = J RU (0). 

As we have mentioned in the introduction, some caution is required before regarding this 
scheme ready for implementation. In practice it will be difficult to achieve perfect degeneracy 
of all states. Thus the question is imposed to which extend incomplete degeneracy of the 
qubit states is permissible. Clearly, the adiabatic condition requires the inverse operation 
time r op to be smaller than the minimum energy difference to the neighboring states: r^ 1 <C 
mm 5 Ec, Jj, Jx- On the other hand, if the degeneracy is not complete and the deviation 
is of the order e one can show by modifying the derivation of Eq. (4) in Ref. |13|| that for 



e <C r op the holonomies can be realized to a sufficient accuracy. This inequality expresses 



the requirement that the operation time be still small enough in order to not resolve small 
level spacings of the order e. 

There is another important constraint on r op . As the degenerate states in Eq. (Q) are 
different from the ground state of the system, r op must not be too large in order to prevent 
inelastic relaxation. The main origin for such relaxation processes is the coupling to a low- 
impedance electromagnetic environment. We can estimate the relaxation rate by Ti n ~ 
E(R env / Rk) where Rk = h/e 2 is the quantum resistance and E is on the order of the 
Josephson energies E ~ Jj,Jx- Thus it is not difficult to satisfy the condition r op < 
T;" 1 experimentally. In fact, it has been found recently that inelastic relaxation times in 
charge qubits can be made quite large and exceed by far the typical dephasing times due to 
background charge fluctuations p3| , p4| , p5| . 

Both charge pumping and the implementation of quantum computing are related to 
coherent manipulations of charge states. Therefore as a readout one can use the scheme 
developed [|1(| to measure charge qubits. No additional difficulty is forecasted at this level. 

Stimulating discussions with G. Falci, G.M. Palma, E. Paladino, M. Rasetti and P. 
Zanardi are gratefully acknowledged. We wish to thank M.-S. Choi for informing us of his 
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FIG. 1. a)Elementary Josephson network for the investigation of non-Abelian geometric phases. 
Note that an asymmetric SQUID loop cannot be switched off completely JH]]. Since Jm = may 
be desirable for quantum computation the SQUID is designed such that this condition can be 
satisfied; b)Pumping cycle for 3 islands. Starting from P = (1/2, 1/2) and adiabatically following 
the drawn path, the gate U = e % ^ (Ty can be achieved. 
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FIG. 2. Geometric realizations of gates U x = ^\ l )u(M ( a ) anc i jj 2 = e^ 2 ^ (b). The structure 
(non-zero elements) of the unitary matrices U\ and U2 is determined by the choice of the plane 
containing the loop and by the starting/ending point of the closed path. Different values of phase 
Si (£2) can be obtained by varying the area enclosed by loops C\ (C2). 
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FIG. 3. Inter-qubit coupling for the implementation of the gate U3 and its geometric realization. 
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